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Abstract 



We present a calculation of diffractive charm jet production at hadron 
colliders in perturbative QCD based on the two-gluon exchange model. Dif- 
fering from the previous calculations, we abandon the use of the effective 
color-singlet vitual gluon simplification, and use the real gluon in the calcu- 
lations for the partonic process. We find the final result is free of the linear 
singularities due to the small transverse momenta of the exchanged two glu- 
ons. In the leading logarithmic approximation (LLA) in QCD, this process 
is related to the off-diagonal gluon density in the proton. As a result, this 
process may provide a wide window for testing the two-gluon exchange model, 
and may be particularly useful in studying the small x physics. In compar- 
ison, the diffractive bottom jet production is also discussed, and is found to 
be important only for large transverse momentum of the jet. 

PACS number(s): 12.40.Nn, 13.85.Ni, 14.40. Gx 
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I. INTRODUCTION 



In recent years, there has been a renaissance of interest in diffractive scattering. These 
diffractive processes are described by the Regge theory in terms of the Pomeron (IP) exchange 
The Pomeron carries quantum numbers of the vacuum, so it is a colorless entity in QCD 
language, which may lead to the "rapidity gap" events in experiments. However, the nature 
of the Pomeron and its interaction with hadrons remain a mystery. For a long time it had 
been understood that the dynamics of the "soft Pomeron" was deeply tied to confinement. 
However, it has been realized now that much can be learned about QCD from the wide 
variety of small- x and hard diffractive processes, which are now under study experimentally. 
Of all these processes, the diffractive heavy quark and quarkonium production have drawn 
specially attention, because their large masses provide a natural scale to guarantee the 
application of perturbative QCD In the framework of perturbative QCD the Pomeron 

is assumed to be represented by a pair of gluons in the color-singlet state. An important 
feature of this perturbative QCD model prediction is that the cross section for the diffractive 
processes is expressed in terms of the off-diagonal gluon distribution in the proton || . 

So far the previous studies are focused on the diffractive processes at ep collider (pho- 
toproduction and DIS processes), we should expect that the two-gluon exchange model can 
also be used to describe the diffractive processes at hadron colliders. Recently, we have 
extended the idea of perturbative QCD description of diffractive processes from ep colliders 
to hadron colliders 0, in which the hadronic diffractive J /if} production is calculated in 
the two-gluon exchange model. In this paper, we will calculate the diffractive charm jet 
production at hadron colliders in this two-gluon exchange model. For a theoretical point of 
view the study of diffractive charm jet production has some advantages compared to J /if} 
production, because it advoids the ambiguities associated with the color-octet production 
matrix elements |7|,|6| and retains the sensitivity to the off-diagonal gluon density in the 
proton. 

However, there exist nonfactorization effects in the hard diffractive processes at hadron 
collisions [[8-11]. These effects mainly focus on the following two sides. One is the so-called 



spectator effects |ITJ. The interaction with the spectator quarks can change the probability of 
the diffractive hadron emerging from collisions intact. This implies that the extra interaction 
with spectators makes it less likely for the diffractive hadron to survive ||12|| . Typically, for 
the diffractive processes at the Tevatron, this survival probability is shown to be about 0.1 
1301. The other one is the "coherent diffractive" processes, in which the whole Pomeron 



is involved in the hard process. These coherent diffractive processes have been proved to 
break the factorization assumption of [l3[. 

As shown in Fig.l, the diffractive charm jet production process pp — > cap calculated in 
this paper belongs to these coherent diffractive processes. The whole Pomeron represented by 
the color-singlet two-gluon system emitted from one hadron interacts with another hadron 
to produce the charm jets. In the leading order perturbative QCD, the partonic process 
gp — > ccp is plotted in Fig.2. In the diffractive final states, there are only charm and 
anticharm jets balanced in the transverse momentum distribution. Experimentally, this 
would provide a strong signal for the coherent diffractive processes at hadron colliders. 

The diffractive production of heavy quark jet at hadron colliders has also been studied 



in Ref. [14 . However, our calculation is quite different from theirs. In their calculation 



2 



of the partonic process gp — > ccp, they used an effective color-singlet virtual gluon (as a 
color-singlet probe in their language) to replace the real gluon for simplification. In our 
calculations, we abandon this simplification, and use the real gluon in the calculations of the 
partonic process. This of course will cause more complicated calculations. However, after 
a lengthy calculation which will be shown in the following we find that our result is free of 
linear singularities, and then guarantees the gauge invariance of QCD. 

Another important difference between and our calculations is the definition of the 



coherent diffraction. Following Ref. ||, we call the process in which the whole Pomeron 
participants in the hard scattering process as the coherent diffractive process. Under this 
definition, all of the nine diagrams of Fig. 2 contribute to the coherent diffractive production 
of heavy quark jet. However, in ]l4j they separate their diagrams into two parts: one part 
(the first four diagrams) contributes to the conventional IS model diffraction; and the other 
part contributes to the so-called coherent diffraction. However, in our calculations, the first 
four diagrams in Fig. 2 alone do not contribute a gauge invariant part of the amplitude, and 
their sum will lead to a linear singularity which is not proper in QCD calculations. So, for a 
full calculation, the nine diagrams must be summed together to contribute to the coherent 
diffraction. 

Another important issue of this process is about the gluon (off-diagonal) distribution in 
the proton at small x. As shown in Ref. [[]], in hadron collisions (such as at the Tevatron), we 
can explore the gluon distribution at small x down to 1CU 6 by studying the hard diffractive 
J/if> production. In the diffractive charm jet production, the involved gluon distribution 
function may take its value at the same order as that for J /if) production. So, detecting 
the diffractive charm jet production at hadron colliders is also attractive for the study of 
the small x physics. Furthermore, experimentally the diffractive charm jet production is 
easily detected because the diffractive final states of c and c pair may have large transverse 
momenta distribution. 

The rest of the paper is organized as follows. In Sec. II, we give the cross section formula 
for the partonic process gp —>■ ccp in the leading logarithmic approximation (LLA) QCD. We 
use the Feynman rule method in the calculations. We also check our method to reproduce the 
formula for the diffractive charm jet production in the photoproduction processes previous 
calculated by others ||. The numerical results are given in Sec. Ill for diffractive charm 
jet production at the Fermilab Tevatron. In Sec. IV, we give some discussions about the 
off-diagonal distribution effects on the process calculated here and the J /if) production of 
H. The conclusion is given in Sec.V. 



II. LLA FORMULA FOR THE PARTONIC PROCESS 

As sketched in Fig.l, the cross section for the diffractive charm jet production at hadron 
colliders {jpp at the Tevatron) can be formulated as, 

da(pp — ► ccp) = J dxida(gp — > ccp)g(xi,Q 2 ), (1) 

where x% is the longitudinal momentum fraction of the antiproton carried by the incident 
gluon. g(xi,Q 2 ) is the gluon density in the antiproton, and Q 2 is the scale of the hard 
process. da(gp — > ccp) is the cross section for the partonic process gp — > ccp. 
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In this section, we first give the formula for the partonic process gp — > ccp in the LLA 
QCD. In the leading order of perturbative QCD, there are nine diagrams shown in Fig.2, 
contributing to the partonic process gp — > ccp. The two-gluon system coupled to the proton 
(antiproton) in Fig.2 is in a color-singlet state, which characterizes the diffractive processes 
in perturbative QCD. The two-gluon exchange model itself may not account for the full 
structure of the Pomeron, but it can approximate the description of diffractive processes in 
terms of perturbative QCD. Due to the positive signature of these diagrams (color-singlet 
exchange) , we know that the real part of the amplitude cancels out in the leading logarithmic 
approximation. To evaluate the imaginary part of the amplitude, we must calculate the 
discontinuity represented by the crosses in each diagram of Fig.2. 

The first four diagrams of Fig.2 are the same as those calculated in the diffractive photo- 
production processes. But, due to the existence of gluon-gluon interaction vertex in QCD, 
in the partonic process gp — > ccp, there are additional five diagrams (Fig.2(5)-(9)). These 
five diagrams are needed for complete calculations in this order of QCD. From the following 
calculations, we can see that these five diagrams are important to cancel out the anormalous 
singularity which rises from the first four diagrams to obtain the correct results. 

To calculate the imaginary part of the amplitude A(gp — > cc+p), we employ the Feynman 
rule method |T5| . For a cross check, we also use this method to calculate the diffractive charm 
jet photoproduction process jp — > ccp. In the leading logarithmic approximation, we can 
reproduce the results obtained previously by using the light-cone wave function method [|J . 



A. Kinematics and Sudakov variables 

In the diffractive process gp — > ccp, the final state only contains the charm quark and 
anti-charm quark. We set kx to be the transverse momentum of the charm quark, m c the 
charm quark mass, and we define rrtj. = m 2 + k\. The invariant mass of the cc pair of the 
diffractive system is set to be M\, and xjp is, 

M x 

xjp = — , 2 
s 

where s is the total cm. energy of the gluon-proton system. 

In our calculations, we express the formulas in terms of the Sudakov variables. That is, 
every four-momenta fcj are decomposed as, 

ki = diq + ^p + k iT , (3) 

where q and p are the momenta of the incident gluon and the proton, q 2 = 0, p 2 = 0, and 
2p ■ q = W 2 = s. oti and /3j are the momentum fractions of q and p respectively, k^ is the 
transverse momentum, which satisfies 

k iT -q = 0, k iT -p = 0. (4) 

All of the Sudakov variables for every momentum are determined by using the on-shell 
conditions of the momenta of the external particles and the crossed lines in the diagram. 

In the following, we calculate the differential cross section da/dt at t = 0. Namely, we 
set the momentum transfer squared of the diffractive process gp — > ccp equal to zero, i.e., 
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u 2 = t = 0. (g + u) is the momentum of the diffractive final state (contains charm and 
anticharm jets), and then 

(q + u) 2 = M 2 x . (5) 

With this equation, and the on-shell conditions of the external proton lines, 

(p-u) 2 =p 2 = J (6) 

we can determine the Sudakov variables associated with u as 

M 2 

a u = 0, Pu = ~f-, 4 = 0. (7) 

In the diffractive region at hadron collisions, we know that M x <C s, i.e., /3 U <C 1. So, in 
the following calculations, we set (3 U to be a small parameter, and take the leading order 
contributions, and neglect the terms proportional to (3 U = ^j L . 

a k and p k are determined from the on-shell conditions of the out-going charm quark and 
anticharm quark, 



(k + q) 2 = m 2 c , 
(u - k) 2 = m\. (8) 

By solving the above equations, we can obtain, 



a k (l + a k ) 



M 2 

(3k = -ot k /3 u = -a k - (9) 

s 

Because a k is of order of 1, the value of m\ is the same order as M\. From the last equation, 
we can see that f3 k is much smaller than a k , i.e., \/3 k \ "C \a k \, because f3 u 1. 

For the loop momentum I, we know that the integral of the amplitude over It receives 
large logarithmic contribution from the region 1/R% l T M\ (Rn is the nucleon 
radius) [Qj. That is to say, the dominant contribution of the integration of It comes from 
the region l\ < M\, so 1% is a small parameter compared with M\ and m T . To calculate 
the integration of the amplitude over l T , we can expand the amplitude in terms of l\ and 
take the leading order contributions. 

The Sudakov variable a\ can be determined from the on-shell condition of the bottom 
cross on the proton line in each diagram of Fig.2, i.e., 

{p-l-uf = 0, (10) 

which results in 

CLl = -4. (11) 

s 

Pi is determined from the on-shell condition of the up cross on the charm quark line 
or the gluon line in each diagram. Unlike other variables calculated above the value of Pi 
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is not the same for these nine diagrams of Fig.2, because there are three different on-shell 
conditions of the up crossed lines, 



(k - I - u) 2 = m 2 c , for Diag.l, 3, 5, (12) 
(k + q + I) 2 = m 2 c , for Diag.2, 4, 6, (13) 
(q + l + u) 2 = 0, for Diag.7, 8, 9. (14) 



These three different relations lead to three different values for fa, 



A= 2(AWT) ~^ , for Diag.l, 3, 5, (15) 

a k s 

2(kr,b) + l> for D . ag _ 
(1 + a k )s 
M 2 - I 2 

A = -, for Diag.7, 8, 9. (17) 

s 

To obtain the above results, we have used the approximations: (3 U <C 1, and \/3k\ <C |afc|. 

To end up the analysis of this subsection, we must note that there are two small param- 
eters in the above calculations of the Sudakov variables, 

(3 U < 1, < 1. (18) 

And these two small parameters are the basic expansion parameters in the following calcu- 
lations of the amplitude for the partonic process gp — > ccp. 



B. Expand the amplitude in terms of \\ 

Using the variables induced in the above, we evaluate the differential cross section formula 
for the partonic process gp — > ccp as, 

da{gp^ccp) dM 2 x d 2 k T da k m 2 

It = 1WM|^ (1 + ak) + Ml ] E 11 ' (W) 

where A is the amplitude of the process gp — > ccp. We know that the real part of the 
amplitude A is zero, and the imaginary part of the amplitude A(gp — > ccp) for each diagram 
of Fig.2 has the following general form, 

ImA = C F {T^ J x Ui(k + q)T ,v 3 {u - k), (20) 

where Cp is the color factor for each diagram, and a is the color index of the incident gluon. 

is some 7 matrices including one propagator. F in the integral represents some other 
factor which is the same for each diagram, 

F=y/j(x>,x";l 2 T ), (21) 
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where 



where the function G(x',x"; k^) is the so-called off-diagonal gluon distribution function ||. 
Here, x' and x" are the momentum fractions of the proton carried by the two gluons. It is 
expected that for small x, there is no large difference between the off-diagonal and the usual 



diagonal gluon densities |16| . So, in the following calculations, we estimate the production 
rate by approximating the off-diagonal gluon density by the usual diagonal gluon density, 
G(x',x";Q 2 ) ~ xg(x,Q 2 ), where x = xp = M\js. Later we will discuss the off-diagonal 
parton distribution function effects for the diffractive charm jet production process calculated 
here and J/tj) production. 

The color factors C F are not the same for the nine diagrams, and they are 



C F = 


2 

9' 


for Diag.l, 4, 


C F = 


1 

~36' 


for Diag.2, 3, 


C F = 


1 

4' 


for Diag.5, 8, 


C F = 


1 

~4' 


for Diag.6, 9, 


C F = 


1 

~2' 


for Diag.7, 



(23) 
respectively. 

As mentioned above, to calculate the leading logarithmic results for the partonic process 
gp — > ccp, the amplitude of each diagram must be expanded in terms of l\. From the 
integral of Eq. (^), we can see that the large logarithmic contribution comes from the 
region l\ -C M\. So, the leading logarithmic contribution to the cross section of the 
partonic process comes from the terms in Y ^ which are proportional to l\. Because we are 
interested in the leading logarithmic results of the cross section, we neglect the higher order 
terms of l\ in F^. 

Futhermore, in the integral of Eq. ( PUD the terms in coming from all diagrams must 
be canceled out by each other. Otherwise, their net sum (order of Ij. ) will lead to a linear 
singularity when we perform the integration over l\. The linear singularity is not proper in 
QCD calculations. So, we first observe the amplitude behavior at the order of l F , i.e., in the 
limit of l\ — > 0. In this limit, the for each diagram has the following result, 



r (l) 


_ r (3) 


r (5) 
A* 


r(2) 


- 


— r(6) 


F(7) 


_ F(8) _ 


_r(9) _ 



Ml 
s 2 

*Pr ( 24 ) 



7 



The second terms in T^> and come from the contributions of Pi terms. This is 

because Pi = — ^f- ^ in the limit l T — > 0. To obtain the above results, we have used the 
following identities, 

u(k + q)(ji+ 4- m c ) = 0, 
(ji— ft + m c )v(u — k) = 0, 
u(k + q)(4+ i)v{u-k) =u[(k'+ 4 - m c ) + ft + m c )]v = 0. (25) 



With the values of the color factor Cp for the diagrams given in Eq. fl23|), we find that their 
sum is zero. That is to say, there is no contribution to the amplitude of the partonic process 
gp — > ccp in this order (1^ takes the order of l T ). As mensioned above, this is what we 
expected. If there are only the first four diagrams (the same as those in photoproduction 
processes), their contributions of the order of l T can not be canceled out by each other, 
because their color factors are not the same for Diag.1,4 and Diag.2,3. (In photoproduction 
processes, the color factors for these four diagrams are the same, so their contribution is 
zero at the order of l T .) However, as mentioned above, QCD gauge invariance requires 
the appearance of other five diagrams at the same order of coupling constant of strong 
interaction. The contributions from Diag.5 and 6 cancel out the linear singularity which 
rises from the first four diagrams. The linear singularity from the last three diagrams, 
Diag.7, 8, and 9, are canceled out by each other. So, the amplitude has no linear singularity 
now. 

From the above analysis, we see that the first four diagrams alone cannot give a gauge 
invariant part of the amplitude for the partonic process gp — *> ccp. 

At the next order expansions of T^, l T , the evaluation is much more complicated. We 
first give the expansion results for the propagators in for every diagrams. To the order 
of l\, the expansions of these propagators are, 



9\ = 7-77T, #4 



a fe M£ {l + a k )M x 

' L 1 ™2 ~ l 1^2 ' 



1 2{k T ,l T ) l\ A{k T ,l T ^ 2 

92 a k M 2 x [ m 2 m 2 (m 2 ) 2 



2(k T ,l T ) l\ , 4(k T ,l T ) 2 



9s = - 7TT-T T772 I 1 + —ZT2 ZZjT + 



yl + a k )M x m\ m\ (ml) 2 

1 r 2(k T ,l T ) , (l + a k )l 2 T , A{k T ,l T f, 



95 ~ M\ [ a k M\ a k M\ (a fc M 



x) 

2(k T ,l T ) a k l 2 T 4(fc T ,/ T ) 2 



96 = TTTi 1 ~ / 1 — T772 + h 7^X772 + 



M 2 x l (l + a k )M 2 x (l + a k )M 2 x ((1 + a k )M x ) 2i ' 

111 

97 = -Tn-. 9s = —9b, 99 = ~— 96- (26) 

M x a k 1 + a k 

Apart from the propagator expansions, the 7 matrices in also contain the l\ terms. 
They are mostly coming from the Sudakov variables (expressed in l T ) of momentum /, i.e., 

oti = — , Pi, and the slasher / T . Furthermore, the 2-dimensional product (fa, It) also 
contributes l\ terms. After integrating the azimuth angle of It, we will get the following 
results, 
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J d 2 l T {k T ,l T f = \\ dl 2 T k T l 2 T , 
J d%(k T , It) h=\J dl 2 T AtI t - (27) 

To test the validity of our procedure, we first evaluate the expansions of T M for the 
first four diagrams, and use these expansions to calculate the amplitude for the diffractive 
photoproduction processes jp — > cap. To do this, we only need do the following changes in 
the amplitude formula Eq. (|20|), 

To the order of l T , the expansions of T„ for the first four diagrams are: 
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r (i) = r (4) = Jr_ 

^-(«fcS i> ^t7m + (1 + ackjsjv ftr fy], 

oikm T 

t— — i tip. - t— — ^ i> - hjw ^ ak ^ 



r (3) 



+ n ~ 1 x 2 (afeg J* #t7/i + (1 + a*)^ (29) 

Using these expansions, together with Eqs. (Effi), (|T^) and (p0[), we obtain the cross section 
formula for the diffractive photoproduction process 7p — > ccp, which is consistent with the 
LLA result calculated via light-cone wave function approach (for Q 2 = 0) [f|. 

For the last five diagrams, the expansions of V ^ are more complicated. To simplify our 
calculations, the contributions of these five diagrams including the color factor Cp are added 
together, and decomposed into several terms as follows. 

The terms (defined as a-term) coming from the slasher 1/p in the 7 matrices, for all these 
five diagrams to sum up together, are 



M-x Aa k (l + a k ) 



i 2 

The terms (6-term) from a\ = —jjt are 

C F Tf = * [(4 + *. J &, fa ~ ^7, 4 i ~ V il, ~ f*7,]. (31) 

4M X a k {l + a k ) a k 1 + a k 2 

The terms (c-term) from f3i are 



CpL » " 2M 2 PM m 2 2 {V ml ] a k (l + a k ) 1 {6 ) 
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The terms (d-term) from propagator expansions are 

r r( rf ) l l ^2, r Qfe ~ C 1 + Qfc) 3 _ ^ + a ^ 

p A{Mlf S ^ [ a k (l + a k ) M\ a 2 (l + a k y J ' ^ 

The last term comes from such terms which are proportional to "(A>r> ^r) /t"- They are 
(e-term) 



C ' Fr ? ) = 77^2^" [~ ~ + ftfc ( 2 «fcg #t + s 7m #t ^- fa i>Vp, - &kS fap^ 
i{M x y a k 



1 + Oik 



[2{l + a fc )s tfk^-s p 1 ^ t7m + ^ # T # M - (1 + a k )s fap„)\. (34) 



Adding up all of the above a, b, c, d, e terms, and those terms coming from the first 
four diagrams in Eq. (|29|) , we get the amplitude squared for the diffractive process gp — > cap 
as, averaged over the spin and color freedoms, 

EW= or/JW 2u m ^ 2 K^ - 8k 2 T m 2 c )(10M 2 x - 27m 2 T ) 2 (xg(x, Q 2 )) 2 , (35) 

y\ 1V1 X) \ m T) 

where the scale Q 2 is set to be Q 2 = m\. The factorization scale in the gluon density is very 
important because we know that the parton distributions at small x change rapidly with this 
scale. Here, we choose the scale to be m\ which is typically used in the calculations of the 
heavy quark production processes. In Ref. [14]], they used different scales for their two parts 



of their diagrams for the partonic process. However, as discussed in the above calculations, 
all of the nine diagrams must be summed together to give a gauge invariant amplitude for 
the coherent diffractive charm jet production at hadron colliders. So, the scales of the nine 



diagrams must be the same. The separation in [14]] is not gauge invariant. 

Finally, the cross section for the partonic process gp — > cc+p, in the LLA of perturbative 
QCD, is 

da(gp —> cc + p) dM\dk\ 7r 2 a^ m 2 

It '' =0 = [M\m\Y 128 x 9 J M 2 X (M 2 X - Am 2 T ) 



(M 2 x m 2 T - Sk^m 2 c )(10M x - 27m 2 T ) 2 (xg(x, Q 2 )) 2 . (36) 
The integral region of k\ satisfies the following constrain 

4 < - m\. (37) 

From Eq. (|36|) , we can see that the cross section of diffractive light quark (high fa) jet 
hadropro duct ion vanishes as m 2 (m q is the quark mass) in the leading logarithmic approxi- 
mation of perturbative QCD (see also 



III. NUMERICAL RESULTS 

With the cross section formula Eq. (^) for the partonic process gp — > ccp, we can 
let the cross section of diffractive production at the hadron level. However, as mentioned 
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above, there exist nonfactorization effects caused by the spectator interactions in the hard 
diffractive processes in hadron collisions. Here, we use a suppression factor Ts to describe 



this nonfactorization effects in the hard diffractive processes at hadron colliders JT0[. At the 



Tevatron, the value of J-'s may be as small as JF S ~ 0.1 [[Hjjn]]. That is to say, the total 
cross section of the diffractive processes at the Tevatron may be reduced down by an order 
of magnitude due to the nonfactorization effects. In the following numerical calculations, 
we adopt this suppression factor value to evaluate the diffractive production rate of charm 
jet at the Fermilab Tevatron. 

Our numerical results are plotted in Fig. 3 to Fig. 7. From Eq. (PE|), we can see that the 
cross section is sensitive to the gluon density in the proton. In fact, the cross section for 
pp — > ccp behaves as 

da^g(x u Q 2 )(g(x,Q 2 )) 2 } (38) 

where x\ is the longitudinal momentum fraction of the proton (or antiproton) carried by the 
incident gluon. So, the cm. energy of the gluon-proton system is s = xiS, where S is the 
total cm. energy of the proton and proton (antiproton) system (e.g., S = (1800GeV) 2 at 
the Tevatron). Then, x = xp = M\j ' s = M^/x\S. We can see that at the Tevatron, x may 
be lowered down to the 10~ 6 ~ 10~ 5 level (varying with the lower bound of the transverse 
momentum kx for the observation). 

In the high energy diffractive processes, we know that the relation M\ <C s = x\S 
must be satisfied (in our calculations, we set M\j 's < 0.1). Together with Eq. (|37|) , these 
constrains will give the integration bound of the arguments xi, M|, and k\. 

In our calculations, the charm quark mass is set to be m c = l.bGeV, and the scale of 
the running coupling constant of strong interaction is set to be Q 2 = m T . For the gluon 
distribution function, we choose the GRV NLO set 



In Fig. 3, we show the M\ dependence of the differential cross section d 2 a/dtdMx\t=o 
for the diffractive charm jet production at the Tevatron. In Fig. 4, we plot the cross section 
da/dt\t=o as a function of lower bound of the transverse momentum kx- From this figure, we 
can see that the cross section is sensitive to the lower bound of kx- It decreases rapidly as 
&Tmm increases. In this figure, we also plot the cross section for bottom quark jet diffractive 
production by using the same formula Eq. (|36| ) except changing the quark mass from m c = 
1.5 GeV to m& = 4.9 GeV. The difference of these two curves shows the dependence of the 
cross section on the quark mass. From Eq. fl3"S|), we can see that at large kx (kx 3> f^b), 
the cross section is proportional to the quark mass, so the cross section of bottom jet 
is larger than that of charm jet. However, due to the existence of the power factor of 
l/m T = l/(k T + m 2 ) in the cross section formula, at small kx smaller quark mass will result 
in larger cross section. So, the cross section of charm jet is larger than that of bottom jet 
at small kxmm- 

In Fig. 5-7, we show the sensitivity of the cross section to the gluon distribution in the 
proton. In Fig. 5, we select xi m i n as argument to show the behavior of the differential cross 
section da/dt\ t =o as a function of xi m i n . We can see that most contribution comes from 
X\ > 10~ 2 region. In Fig.6, we plot the cross section as a function of a; max , which shows that 
most contribution comes from x ~ 10 -3 — 10 -2 region. In these two figures, for both we set 
the lower cut of the transverse momentum kx to be kx m m — 5 GeV. From Fig.4, we see that 
the lower bound of transverse momentum kxmm largely affects the total cross section of this 
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process. Furthermore, it is expected that /cTmm can also affect the relative importance of 
different x region contribution to the total cross section. This effect can be seen from Fig. 7. 
In this figure, we plot the ratio o /a tot as a function of x max , where a to t is the total cross 
section after integrating over all x region, while a is the cross section with the upper bound 
s max for x. From left to right, these three curves correspond to three different lower bound 
of transverse momentum, /cTmin = 0, 5, 10 GeV respectively. From Fig. 7, we see that for 
smaller value of fcTmin, the dominant contribution comes from smaller x region. 

As stated in the section of Introduction, the nine diagrams in our calculations all con- 
tribute to the coherent diffractive charm jet production at hadron colliders. So, the above 
calculated numerical results from Figs. 3-7 are all for the coherent diffractive production. 



IV. THE OFF-DIAGONAL GLUON DISTRIBUTION EFFECTS 



The off-diagonal parton distribution function is a hot topic in recent years ||19|| . Typically, 
the off-diagonal distribution has two different regions. One is x' > xjp (x" > 0, because 
x' — x" = xjp), and the other is x' < xp. For the off-diagonal gluon distribution function 
G(x' , x"; Q 2 ), x' and x" are the momentum fractions of the proton carried by the outgoing 
gluon and the returning gluon respectively. In the first region, because x" > 0, G(x', x"; Q 2 ) 
is similar to the usual gluon distribution function g(x;Q 2 ). However, in the second region 
(x" < 0), G(x',x";Q 2 ) is rather a distribution amplitude of taking two gluons out of the 



proton, which is similar to a meson's wave function. Previous studies |L6[ have shown 
that in the first region, the off-diagonal gluon distribution function is not very different 
from the usual diagonal gluon distribution at small x. So, in this region we can safely 
approximate the off-diagonal gluon distribution function by the diagonal distribution. On 
the other hand, in the region of x' < Xjp, where the off-diagonal gluon distribution function 
is distinctively different from the traditional diagonal gluon distribution function P0| , the 
off-diagonal distribution function effects may be important. 

From Fig. 2, we find that there are three different values for x' (x") because there are 
three different values for j3i, 

x' = (3i + (3 U = xjp + 2 ( kT ^- l T for Diag-1> 3; 5) 

a k s 



x F + 2 ^ + 1 t for Diag . 2; 4 6j 
(1 + a k )s 



for Diag.7, 8, 9. (39) 



When integrating the amplitude over It, x' may be smaller than xjp in some region, where 
the off-diagonal gluon diatribution function G(x f , x"; l T ) will take its value in the second re- 
gion. Especially, for Diag.7-9, the dominant contribution to the integration of the amplitude 
comes from the region where x' is mostly smaller than xjp, because the large logarithism 
contribution of the integration comes from the region of x' xjp (l T M^). So, we can 
expect that the off-diagonal effects on this process may be important. 

In our previous calculations of the diffractive J / ip production at hadron colliders, we also 
neglect the off-diagonal gluon distribution function effects. However, from Fig.2 of Ref. ||, 
we see that x' has two different values, 
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for Diag. d, e. (40) 

When integrating the amplitude over l^, the off-diagonal gluon distribution functions for 
Diag. a, b, c take their values in the region of x' < xjp, in which we can safely use the usual 
gluon distribution instead. However, for Diag.ci, e, the off-diagonal gluon distributions dom- 
inantly take their values in the region of x' < xjp (because the large logarithmic contribution 
to the amplitude comes from the integral region l\ <C Ml). So, the off-diagonal distribution 
effects in the diffractive J/ip production at hadron colliders are also expected to be very 
important. 

To evaluate the off-diagonal distribution effects on the diffractive processes (charm jet 
and J /if} production at hadron colliders) discussed above, we must integrate the amplitude 
including the differential off-diagonal gluon distribution function f(x',x";l^) in Eq. ( |22|) 
(with the evolution scale at Q 2 = 1?) over l\ numerically. This task is difficult and beyond 
the scope of this paper. Further investigation in this direction is in progress. 

V. CONCLUSIONS 

In this paper, we have calculated the diffractive charm jet production at hadron colliders 
in perturbative QCD by using the two-gluon exchange model. Differing significantly from 
other models, our coherent diffractive amplitude for the partonic process is gauge invariant 
and contains no linear singularities as l\ — > 0. We give the formula for the cross section of 
this process in the leading logarithmic approximation. By neglecting the off-diagonal gluon 
distribution effects, we have estimated the production rate of diffractive charm jet production 
at the Tevatron. The comparison of production rates between the charm jet and the bottom 
jet has shown the dependence of the cross section on the quark mass, which is a distinctive 
feature of these processes. We have also discussed the off-diagonal gluon distribution effects 
on this process and the hard diffractive J/ip production previously calculated ||. 

In conclusion, we have calculated the heavy quark jet production in the diffractive pro- 
cesses, which belongs to the coherent diffractive processes at hadron colliders. Detecting the 
heavy quark jet in hadron collisions would provide strong signals for these coherent diffrac- 
tive processes. We note that the process calculated in this paper is gluon initiated process, 
which is sensitive to small- x gluon distribution in the proton. Some other processes, such 
as W ± , Drell-Yan processes, are quark initiated processes, which would be sensitive to the 
large-x quark distribution function in the proton ||. Work in this direction is in progress. 
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Figure Captions 



Fig.l. Sketch diagram for the diffractive charm jet production at hadron colliders in pertur- 
bative QCD. 

Fig. 2. The lowest order perturbative QCD diagrams for partonic process gp — > ccp. 

Fig. 3. The differential cross section d 2 a/dtdMx\ t=0 at the Fermilab Tevatron as a function 

of M|-. 

Fig. 4. The differential cross section da/dt\ t=0 after integrating over k\ (for k\ > k^ min ) as a 
function of k Tniin , where m c = 1.5 GeV for the charm jet and m& = 4.9 GeV for the bottom 
jet. 

Fig. 5. The differential cross section da/dt\t=o as a function of xi m j n , where xi m i n is the lower 
bound of x\ in the integration of the cross section. 

Fig. 6. The differential cross section da/dt\ t =o as a function of a; max , where x max is the upper 
bound in the integration of the cross section. 

Fig. 7. The ration a/a tot as a function of x max , where a tot is the total cross section after 
integrating over all x region. From left to right, the three curves correspond to three lower 
bounds of transverse momentum, /cTmin = 0, 5, 10 GeV respectively. 
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